A critical note on one-eighth law in 2D turbulence 
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We critically revisit the various attempts to prove one eighth law in 2D turbulence and reconcile 
them. Herein, the one-eight law has been proved for unforced 2D incompressible high Reynolds 
number turbulence. An exact expression of the time derivative of two-point second order velocity 
correlation function is also derived for the enstrophy cascade dominated regime. 
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The phenomenological model of Kolmogorov[l] remains the cornerstone of our knowledge about the theory of 
3D incompressible high Reynolds number turbulence. The most important non-trivial result that it yielded is the 
four-fifths law: 




[v(f + I) — v(r)\ .— 
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where v is the velocity field and e is the rate per unit mass at which energy is being transferred through the inertial 
^ . range. This rare exact result has stood the rigorous tests put forward by experiments and numerics. Theoretically this 
result could be derived in two ways [2]: i) the original Kolmogorov method explained in details in the fluid dynamics 
text due to Landau and Lifshitz[3]. (There is no external forcing in this approach and the equality of dissipation rate 
^-y and forcing rate for the energy is never enforced), and ii) a field-theoretic method. (Here an external forcing maintains 
\ a steady state turbulence) . The field theoretic technique invokes the so-called "dissipation anomaly" for deriving the 
, Ph - law. 

Even in 2D turbulence one would like to seek such a classic law concerning two-point third order structure function 
t-H ■ S3 using one or the other of the above-mentioned methods. Recent attempts to do so by Lindborg[4], Bernard [5] 
1 and Chakraborty[6] amounted to some revealing results. The equivalence between the two techniques has also got 
highlighted by virtue of these results. However, one point of confusion has remained. It concerns the inability 
to derive one-eighth law in the enstrophy cascade dominated regime in unforced 2D turbulence using Kolmogorov- 
Landau method [6]. It may be contrasted with the fact that the one-eight law seems to be correctly derived in 2D 
. ■ forced turbulence[4, 5]. 

Frisch[7] has observed that as certain mathematical steps leading to the proof of one-eighth law in forced 2D turbulence 
are of purely kinematic nature and consequently they should also lead to the one-eighth law in unforced version of 2D 
turbulence. If such is the case then it would modify certain conclusions of reference[6] and would extend the validity 
of one-eighth law into the regime of unforced 2D turbulence. We shall see shortly in what follows that this indeed is 
the case. 

It is a well-established fact that there exists a direct-cascade of enstrophy in two-dimensional (2D) turbulence. One 
defines total enstrophy as T = i J all spacc u> 2 d 2 p where uj = d x v y — d y v x is the vorticity in the Cartesian coordinates; 

v being the velocity field. As we shall consider incompressible fluids only (V.v — 0), we shall take the density to be 
unity and let p take over the task of representing position vector in 2D plane. The enstrophy flows through the inertial 
range and gets dissipated near the dissipation scale. Using the antisymmetric symbol e a p that has four components, 
viz. e\\ = £22 = and £ 12 = —£21 = lj one may define the mean rate of dissipation of enstrophy per unit mass as: 

r) = v(Voj.Vu)) 

=*> 7] = ueraeeiiiidrd^a^ded^p)) (1) 

Here, angular brackets denote an averaging procedure which averages over all possible positions of points 1 and 2 at a 
given instant of time and a given separation. Now, if v\ and v% represent the fluid velocities at the two neighbouring 
points at pi and p 2 respectively, one may define rank two correlation tensor: 

B a p = {{Via - «la)(«2/3 ~ Vip)) ( 2 ) 
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For simplicity, we shall take a rather idealised situation of turbulence flow which is homogeneous and isotropic on 
every scale, a case achievable in practice in a vigorously-shaken-fluid left to itself. The component of the correlation 
tensor will obviously, then, be dependent on time, a fact which won't be shown explicitly in what follows. As the 
features of local turbulence is independent of averaged flow, the result derived below is applicable also to the local 
turbulence at distance p much smaller than the fundamental scale. Isotropy and homogeneity suggests the following 
general form for B a p 

B afi = A l (p)5 afi + A 2 (p)p a p° p (3) 

where A\ and A2 are functions of time and p. The Greek subscripts can take two values p and _L which respectively 
mean the component along the radial vector p and the component in the transverse direction. Einstein's summation 
convention will be used extensively. Also, 

P = P2~Pl, p° a = p a /\p\, P° p = l, P°L = 

using which in the relation (3), one gets: 

B a p = B ±± (6 a p - p° aP %) + B pp p° a p% (4) 

One may break the relation (2) as 

B a p = (viaVip) + (v2aV2f]} ~ (vi a V2p) ~ {v 2a V\p) (5) 

and defining b a p = {v\ a V2p) one may proceed, keeping in mind the isotropy and the homogeneity, to write 

B a p = (v 2 )S a f3 - 26 Q/3 (6) 

Let's again define: 6 a /3 i7 = {vi a v 1 pv 2l )- Invoking homogeneity and isotropy once again along with the symmetry in 
the first pair of indices, one may write the most general form of the third rank Cartesian tensor for this case as 

b afin = C(p)S aP p , ( + D(p)(S lP p a + S ai p$) + F(p)p a p$ P ; (7) 

where, C, D and F are functions of p. Incompressibility dictates: 

d d 

a b a p.y = q Oafj.j = (8) 

op2 7 ap 1 

„ /r C r 2D 2D' 2D F' F 

=4- C a p H S a fj H 5 a p H p a pp ^PaPfi + ~^PaPP + ~^PaPP = (9) 

p p p l p A p z p 6 

Here prime (') denotes derivative w.r.t. p. Putting a = [5 in equation (9) one gets: 

^ „ constant , N 

2C + 2D + F= = (10) 

P 

where, it as been imposed that b a p^ should remain finite for p = 0. Again, using equation (9), putting a ^ (3 and 
manipulating a bit one gets: 

D = -\(pC' + C) (11) 

using which in relation (10), one arrives at the following expression for F: 

F = pC -C (12) 

Defining 

B a pj = ((v 2a - v la )(v 2 p - vip)(v 27 - «i 7 )} (13) 
= 2(b a/3n + b 7/3 . a + b a7i p) (14) 

and putting relations (11) and (12) in equation (14) and using relation (7), one gets: 

B aPl = -2 P C'(5 a pp° + 8 lP p° a + 6 ai p°p) + 6(pC - C)p a p p p° y (15) 
=> S 3 = B ppp = -6C (16) 
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which along with relations (11), (12) and (7) yields the following expression: 

Kp n = + ^(pS' 3 + Ss)(5 10P ° a + 5 a ,p^) l -{pS' 3 - S 3 )p° aP ° p° (17) 

Navier-Stokes equation gives: 

d 

-Q t v ia = -viydi^a - diaPi + vd l7 d l7 v la (18) 

-Q t v W = - V 2 1 d 2l v 2 p - d 2 pp 2 + vd 2l d 2l v 2 f3 (19) 
Multiplying equations (18) and (19) with v 2 p and v\ a respectively and adding subsequently, one gets the following: 

d 

-Q- t {vi a v 2 p) = -d l7 {v l7 v la v 2 p) - d 2l {v 27 v la v 2 p) 

-di a (piv 2f3 ) - d 2 p{p 2 vi a ) + vdi 7 di 7 (vi a v 2 i3) 

+vd 2l d 2l (vi a v 2 p) (20) 

Due to isotropy, the correlation function for the pressure and velocity ((p\v 2 )) should have the form f(p)p/\p\. But 
since d a (piv 2a ) = due to solenoidal velocity field, f(p)p/\p\ must have the form constant x (p/\p\ 2 ) that in turn 
must vanish to keep correlation functions finite even at p = 0. Thus, equation (20) can be written as: 

d 

-Q-b a f3 = 9 7 (6 Q7i/3 + b^.a) + 2vd 1 d 1 b aP (21) 

For isotropic and homogeneous turbulence, the condition of incompressibility gives the easily derivable well-known 
result: 

4:d^b a j a — d~yB aa ry (22) 
Defining W = {uo\uj 2 ) and noting that W = —dsdsb aa , we get from relations (21) and (22): 

dW 1 



dt 2 



d 5 d 5 {d 1 B aai )-2ud 5 d s W (23) 



Again, if one defines £1 = ((ui 2 — u\)(u) 2 — t^i)), for homogeneous isotropic turbulence one may write O = 2(uj 2 ) — 2W . 
So, equation (23) can be manipulated into the following: 

1 f)Q f)(bJ 2 \ 1 

2~dt~ dt = 2 dsd5d ~< Ba ^ + vdsd&n ~ 2vdsds ^ 2 ) ( 24 ) 

=> dsd s d 7 B aai = 4?7 (25) 
1 

V 



B aa - ( = 7 VP 3 (26) 



Here, we have assumed ^ to be relatively negligible and let v — ► so that the terms proportional to v vanish. Also, 
we have recalled that ^ \ = _^ From equations (7) and (14), and the condition of incompressibility, it readily 
follows that B±± p = £j^B ppp putting which in expression (26) and integrating subsequently (keeping in mind that 
Bppp shouldn't blow up at p = 0), we arrive at: 

B PPP = +\vP 3 (27) 

This is the one-eighth law for the unforced 2D incompressible turbulence proved using the Kolmogorov-Landau 
approach. 

Let us go back to equation (21). Using expressions (6) and (17), one can rewrite the equation as: 

19,,, 1 3 „ „ „ , 19,.. \ v d ( dB 00 \ , . 

2d& > - 28^ = ^> > + W*BP {p BpPp) " 'POP- V^V) (28) 

As we are interested in the enstrophy cascade, the first term in the R.H.S. is zero due to homogeneity and the first 
term in the L.H.S. is zero because of energy remains conserved in 2D turbulence in the inviscid limit (and it is the 
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high Reynolds number regime that we are interested in); it cannot be dissipated at smaller scales. Also, as we are 
interested in the forward cascade which is dominated by enstrophy cascade, on the dimensional grounds in the incrtial 
region B pp (as it may depend only on 77 and p) may be written as: 



d_ 

dt 



-B pp = A VP 2 (29) 



where A is a numerical proportionality constant. Hence, using the relation (29), the equation (28) reduces to the 
following differential equation: 

1 d , o N v d ( dB OD \ A 9 , . 



which when solved using relation(16) in the limit of infinite Reynolds number {v — > 0), one gets 

B PPP = -^ P 3 (31) 



Comparing it with the expression (27) for the two-point third order velocity correlation function for the isotropic and 
homogeneous 2D decaying turbulence in the inertial range of the forward cascade, we determine the value of A to be 
-1/4. Therefore, relation (29) yields 

-B pp = - iW 2 (32) 

This, to the best of our knowledge, is yet another exact new result that has to be verified experimentally and 
numerically to test its validity. 

Having answered the confusion mentioned earlier regarding the one-eighth law, we turn critical to yet another issue 
that was not addressed in any of the references mentioned in the beginning. We have been careless enough to assume 
the existence of r\ when v — > 0. If 77 doesn't exist, the one-eighth law is yet again in jeopardy. It is really unfortunate 
for the law that it has been rigourously proved [8, 9] that enstrophy dissipation is not possible for any 2D Eulcr 
solutions with finite enstrophy. Thus, rj may exist in the inviscid limit only when one takes rather ill-defined initial 
conditions for which the total initial enstrophy is infinite. In view of this, one must not be surprised at all if numerics 
and experiments fail to uphold the one-eighth law in many a situation. The law of 2D turbulence, therefore, doesn't 
enjoy the same classic status as the Kolmogorov law of 3D turbulence. 
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